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Non-Abelian color fields from relativistic color charge configurations 

in the classical limit 
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We study the dynamics of color fields as generated by simple configurations of relativistic particles 
with Abelian and non-Abelian (5(7(2)) charges in the classical limit. We find that chromodynamic 
(non-Abelian) systems generally show Coulomb-like features by analogy with electrodynamics. A 
peculiar feature in the non-Abelian case is the additional strength of the chromoelectric and chromo- 
magnetic fields caused by the contribution of changing the color charge. This change of color SU{2) 
charges results in a rotation of the color vector which is getting very fast at close partonic distances. 
The presence of this non-Abelian additional term in the chromoelectric and chromomagnetic fields 
creates also a 'color charge glow', a significant contribution to the interaction of the partons due 
to the time dependence of the color charge. This situation is especially relevant to the very initial 
phase of ultra-relativistic heavy-ion collisions, where the initial partonic state is governed by high 
(non-equilibrium) parton densities and strong local color fluctuations. 

PACS numbers: 25.75.-q, 24.85. +p, 12.38.Mh 



I. INTRODUCTION 



Ultra-relativistic nucleus-nucleus collisions at Rela- 
tivistic Heavy-Ion Collider (RHIC) or Large-Hadron- 
Collider (LHC) energies provide an opportunity to ex- 
plore strongly interacting QCD matter in- and out-off 
equilibrium. The experiments at the RHIC and the LHC 
have demonstrated that a stage of partonic matter is pro- 
duced in these reactions which is in approximate equi- 
librium for a couple of fm/c [ll-Q. Contrary to early 
expectations - guided by perturbative QCD (pQCD) cal- 
culations - the medium showed a collective response of 
the strongly interacting plasma (sQGP) with interaction 
rates beyond those of high density hadronic matter . 
In fact, viscous hydrodynamical calculations are success- 
ful in describing the collective response of the partonic 
medium which is dominantly reflected in the differen- 
tial azimuthal angular distribution of the final (observ- 
able) hadrons |9l-[ll|. These angular distributions can 
be characterized by the Fourier coefficients w„ that spec- 
ify the strength of the n-th harmonics in the angular 
distributions. The observation of sizeable uneven coef- 
ficients V3,V5 showed the presence of initial-state fiuc- 
tuations which are due to random positions of the nu- 
cleons in the targets and random inelastic hard interac- 
tions in the primary phase of the nucleus-nucleus colli- 
sions. This comes about as follows: the passage time of 
even central Au-|-Au (or Pb+Pb) at RHIC energies is no 
longer than 0.2 fm/c at top RHIC energies and below 
0.02 fm/c at present LHC energies. During these short 
times the wave- functions of the nucleons in the target are 
frozen and even the partonic wave-functions of the indi- 



vidual nucleons are approximately frozen so that configu- 
rations with high fluctuations in the energy density show 
up in individual events (l^ - [l6j . In fact, when account- 
ing for such initial state fluctuations all harmonics w„ 
in the azimuthal angular distributions may be described 
in event-by-event calculations within the hydrodynamic 
framework fl7{ or within transport approaches like the 
Parton-Hadron-String Dynamics (PHSD) theory (lH - 

The partonic medium (close to equilibrium) is char- 
acterized by a very low ratio of the shear viscosity to 
the entropy density rj/s ~ 0.1-0.2 which is close to the 
lower bound of 1/ (47r) [l^ . It is presently accepted that 
rj/s has a minimum as a function of temperature close to 
Tc ~ 160 MeV and increases fast in the hadronic phase, 
i.e. with decreasing temperature [20l - [2^ . Furthermore, 
the ratio of the bulk viscosity to entropy density is ex- 
pected to have a maximum close to Tc |23l425l |. 

However, the (dissipative) hydrodynamic calculations 
start at delayed times of 0.2 - 1.0 fm/c from the initial 
impact of the heavy ions since approximate equilibrium 
has to be achieved by presently unknown mechanisms 
[26j . Initial fluctuations in the hydro calculations have 
to be imposed by independent Glauber model simula- 
tions or Color Glass Condensate (CGC) initial conflgura- 
tions [131 • Initial conditions very similar to the Glauber 
model are included by default in the PHSD transport 
approach that has been tested successfully so far for a 
variety of observables in nucleus-nucleus collisions from 
lower Super-Proton-Synchrotron (SPS) up to LHC en- 
ergies [28-31], includingelectromagnetic probes such as 
e^e~ or pairs [32|. The event-by-event PHSD 
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energies show that apart from local fluctuations in the 
energy density also sizeable fluctuations in "color" oc- 
cur although the total system is always color neutral. 
Furthermore, the formation of color dipoles is observed 
in the microscopic transport calculations, which implies 
that color forces should play an important role at least 
during the early non-equilibrium stage (cf. Ref. [33|). 
Such color forces occur for all centralities of the collision 
and might be even stronger than the average forces in 
the partonic medium reflected in the ratio rj/s since the 
shear viscosity is inversely proportional to the average 
transverse cross section between the constituents. 

Our present study is, furthermore, motivated by re- 
cent investigations of peripheral ultra-relativistic heavy- 
ion collisions at the RHIC and the LHC, for which a large 
electric charge eZ of the colliding nuclei leads to the gen- 
eration of intense electric and magnetic fields during the 
passage time of the charged 'spectators', as discussed in 
Ref. [3J|. It is speculated that a very strong electro- 
magnetic field - of short duration essentially in the pre- 
equilibrium phase - might have an important impact on 
particle production [35j, [SGj] . Fortunately, a large theoret- 
ical (phenomenological) activity combined with enthusi- 
asm in studying the chiral magnetic effect (CME) and, in 
general, strongly interacting matter in magnetic fields in 
recent years [37] is balanced by a careful analysis of the 
corresponding RHIC and LHC experimental data j!!] on 
the charge separation dependence of azimuthal correla- 
tions and, in particular, correlations with respect to the 
reaction plane. Actually, the measurements of the AL- 
ICE Collaboration [s^ of the pion angular correlations 
are consistent with the qualitative expectations from the 
chiral magnetic effect. However, these results are also 
consistent with local charge conservation in combination 
with a large elliptic flow V2- 

The future ALICE measurements are promising to 
provide an answer to these questions and, possibly, 
to strengthen the standard approach to the satura- 
tion physics [1^ treating the gluon fields and corre- 
sponding scattering cross sections classically (rapidity- 
independent). Such classical field dynamics call for a 
transport formulation and the development of an ex- 
tended transport code for ultra-relativistic collision pro- 
cesses, including the color dynamics of the gluon fields 
(beyond the processes implemented e.g. in PHSD). 
One might believe that the problem of determining the 
proper non-Abelian fields could easily be resolved by av- 
eraging over the color sources and then applying the 
linear BFKL approach (including quantum correc- 
tions) or the non-linear evolution equations by Balitsky- 
Kovchegov [i^ and [i^. However, present experience 
tells us that this way is (at least) inefficient. 

In order to proceed to a solution of this task, we 
are going here to give practical estimates for the space- 
time dependence of classical non-Abelian fields which 
(in line with PHSD) are generated by colored quarks 
at the initial stage of the heavy-ion collision. In this 
respect we will develop an instrumental approximation 



for configurations of two color charges, a color charge 
and a color dipole and two color dipoles which are mov- 
ing along a straight line towards each other. For our 
initial study of non-Abelian fields we simplify the prob- 
lem to the SU(2)-color group. We recall that such color 
configurations are just the elementary basic examples in 
ultra-relativistic heavy-ion collisions in the existing phe- 
nomenological models [28, 44]. The classical character of 
the non-Abelian field means, as usual, that the field op- 
erators are replaced by their average values for the quan- 
tum state, and the off-diagonal elements are neglected. 
Quarks are treated as classical point-like massive par- 
ticles possessing a color charge. We will estimate the 
characteristic values of the classical gluon field strength 
but are not interested in the parton (gluon) distribution 
functions themselves. More precisely, we will specify the 
features characteristic of the non-Abelian nature of fields 
in comparison to solutions for Abelian fields for the same 
coupling strength. 

The actual layout of our study is as follows: In Section 
II, we first recall results from classical (Abelian) electro- 
dynamics considering the retarded electromagnetic field 
created by a moving source. In Sect. HI the particular- 
ities of the non-Abelian color field produced by a color 
source are discussed on the basis of an approximate so- 
lution of the classical Yang-Mills equations. The appli- 
cation of this technique to particular cases of two color 
charges, a charge-color dipole system and two dipoles 
moving in the opposite direction is presented in Sections 
VI, V, VI, respectively. Our conclusions are given in 
Section VII. Some more technical details for the color in- 
teractions and the construction of approximate solutions 
of the classical Yang-Mills equations are shifted to the 
Appendix. 



II. THE FIELD OF A RELATIVISTIC ABELIAN 
CHARGE 

In this Section, we recall results from classical electro- 
dynamics: The field of a point-like charge propagating 
along the trajectory r{t) is described by the (retarded) 
Lienard-Wiechert potential at the observation point 



_ 1 



R~vR 



An 



e V 



R-vR 



(1) 



Throughout the paper we stick to the standard system 
of units with speed of light c = 1, with dimensionless 
electrodynamic e and non-Abelian g charges, i.e., have 
the proper factors of he to keep the proper dimensions of 
magnitudes. Then characteristic distances in the prob- 
lem are of an order of one Fermi and the potentials and 
field strengths are measured in units of m,r and mj, re- 
spectively, where to^ is the 7r-meson rest mass. In Eq. 
([1]), is the O'th component of the potential and A is 
its vector component, v is the particle velocity at some 
retarded time t' which is determined by the distance be- 
tween the observation point and the particle R = \R\ 
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and R = ro — r(t') is the radius-vector from the charge 
position to the observation point r^. It is of importance 
to note that the relation between the retarded and labo- 
ratory time t is 



R = t-t' , 



(2) 



(although notation R' would be more logical). The elec- 
tric and magnetic fields then are given by 

dA 

E = -—-Vip, H = VxA. (3) 



Using Eq. ^ we get 
OA dA dt' 



dt dt' dt 
dt' 1 



dt 1 — vn 
V{vR) = V- 
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with the unit vector n = R/R while v = dv/dt' is 
the particle acceleration retarded in time (see Ref. (45|). 
Combining these relations and using the potentials from 
Eq.(IT]) we get the following results for the electric and 
magnetic fields: 
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In fact the formulae (convenient for practical purposes) 
express the electric and magnetic fields as a function of 
the particle trajectory (at v ^ 0) at the current time 
r — ro — r{t) in the form 
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e r(l — v'^) 



H = vx E . (6) 



An [(r-i,)2 + r2(l--i;2)]3/2 
We also recall that the potentials can then be written as 



An [(ri;)2 + r2(l--y2)]i/2' 
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An [(ri;)2 + r2(l-i;2)]i/2- 



III. THE FIELD OF RELATIVISTIC SU{2) 
CHARGES (CLASSICAL YANG-MILLS FIELDS). 

In order to keep our approach transparent we consider 
here the SU(2) non-Abelian fields only. In the classical 
approximation they are c-number functions which are the 



solutions of the classical Yang-Mills equations. We ana- 
lyze the QCD Lagrangian of the form 



1 



(8) 



where the color vector = (A^,A^,A^) represents 
a triplet of the Yang-Mills fields of different colors 
("isospin"), is the current density of external color 
sources and G^^ = — d^A/^ + gA^ x A^ is the 

gluon field tensor with the covariant derivative acting as 
Dt^f = d^ f + gA^^ X f. The product sign x corresponds 
to the vector product in the color space. The classical 
equations of motion, as known, then read 



(9) 



The most significant difference between these equations 
and the electrodynamic equations is the compatibility 
conditions of the system of Eqs. ^ 



(10) 



that generally implies that the color vector charge is not 
conserved in a sense similar to electrodynamics. Then 
as a suitable solution of the Yang- Mills equations for a 
single particle with constant color charge C one may take 
the potentials of the form 
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Similarly to the electromagnetic case discussed above in 
Section II R = — r{t') is the radius-vector directed 
from the charge position to the observation point Tq. In 
electrodynamics one can easily verify the uniqueness of 
the solutions ^ for a point-like charge (taking into ac- 
count the well-known problem of the singularity associ- 
ated with the self-interaction [46]) because the superpo- 
sition principle for solutions remains valid. Furthermore, 
the general form of the retarded solution is well known. 
In the non-Abelian case we cannot benefit from an anal- 
ysis of the general Yang-Mills solutions but one can con- 
struct approximate solutions with properties similar to 
the solutions of electrodynamics. As mentioned above, 
here we focus on the analysis of some simple examples im- 
portant for practical applications treating the relativistic 
color objects propagating along a straight line. In view 
of the high parton densities in the initial stage of ultra- 
relativistic nucleus-nucleus collisions the average distance 
between the color charges is much smaller than 1 fm such 
that long-range color forces are screened and of minor im- 
portance. It is expected that color Coulomb interactions 
should be dominant in this initial phase. 

As an approximate solution, we consider a superpo- 
sition of the Lienard-Wiechert potentials (|lip in which 
the vector of the particle color charge can change in time 
and should be taken at the retarded time which results 
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exactly from the general form of the retarded solution as 
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(12) 
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As an example, we consider a simple approximation of 
the covariant 4-current with the delta function for a 
point-like particle 

j^(r', t') = [C S{r' - rit)),C v{t) 5{r' - r(0)] , (13) 

where v{t) = r{t) is the particle velocity. In this context 
the compatibility condition ((TO)) for a point-like charge 
becomes 



C = g ^{t,r)-v A{t,r) xC 



(14) 



From these equations it is easy to see that the modu- 
lus of the color charge vector remains constant. In the 
following we take the color charge as a unit vector and 
specify the magnitude of the charge to be given by the 
coupling constant g {ag — .g'^/(47r) = 0.3), thus taking 
the coupling strength out as an independent factor. This 
approximation will be justified below in more detail. 

The consideration of the non-Abelian charge as a func- 
tion of time leads to the generation of additional terms in 
the expressions for the chromoelectric and chromomag- 
netic fields of color point-like charges which stem from 
the derivative of the color charge vector w.r.t. the re- 
tarded time t', i.e. 
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V X Cv = Vt' X Dv 
The results for the non-Abelian chromo-fields are 
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These additional terms look like radiation terms, how- 
ever, one cannot state that exactly because the change 
of the color charge vector does not result from a par- 
ticle movement and such a conclusion needs a further 
special analysis. It is difficult also to compare directly 
these results with the Abelian case since the chromo- 
electric and chromomagnetic field strengths are not ob- 
servable quantities. Accordingly, it is more consistent to 



consider the energy density of chromoelectric and chro- 
momagnetic fields /2 or the Wong force which 
is analogous to the Lorentz force in electrodynamics [47 1 
and serves as a measure of the momentum change of a 
color particle in an external color field. Such a role could 
also be played by a set of covariant functions which de- 
fine the strengths of chromoelectric and chromomagnetic 
fields in the different space-time points (in principle, one 
could determine the covariant gluon-parton distribution 
functions in the following Weizsacker- Williams ansatz). 

In order to get well-settled approximate solutions, we 
need further information that, e.g., can be obtained 
by solving the problem of two slowly moving color 
charges O . Assuming one particle to carry a color 
charge P and the other particle a charge Q we are al- 
lowed to simplify the consistency equations ([TU)) for the 
point-like sources to the following form : 



P = g ip{t,xi) X P , 
Q ^ g (p{t, X2) X Q , 



(16) 



where the dot denotes the time derivative, Xi and X2 
are the positions of the point-like charges at time t and 
tp{t, Xi) is the value of the scalar potential at the charge 
location. An immediate consequence of Eqs. ([T^ is the 
conservation of the modulus of the particle color charge 

P'^ — and Q'^ — 0. An analysis of the perturbation 
series in the coupling constant g and in particle velocity 
v/c shows that if the scalar potential ip is spanned by the 
two color vectors defined above 



tp ^ (fil P + ip2 Q , 



(17) 



(in particular, the first term of the iteration series looks 
similarly) , then the vector potential A (generated by the 
charges) is spanned by the vector product of charges only 



A = a PxQ 



(18) 



In fact, Eqs. (IT71) and (IT51) fix the Coulomb gauge. More- 
over, as a consequence of the vector algebra rules, contri- 
butions to the scalar potential are given only through the 
terms spanned by two charge vectors P, Q. Finally, the 
equation is factorized and we have the following system 
of equations for the potential components ipi, ip2 and the 
vector field a 



DD $ = (5 , 

VxVxa-gj = 0, j = $ JZ)$ 



(19) 



where the column $ 



ip— ^ with ip^ 



V = II <^i,</?2 II, (/5i= Lpx(x2) and (^*2= P2(xx) is the 
potential at the points of the charge location, = 
\8(x - xx),5{x - 0:2)11, Dki = V4/ +gaCki, 
k,l — 1,2 is the covariant derivative and C, J are the 
following matrices : 
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-1.5 -1 -0.5 0.5 1 1.5 2 2.5 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 

z [fm] z [fm] 



Figure 1: (Color online) Schematic view of the lines of force of the vector field a (left panel) and the isolines of the chromomag- 
netic field component Hip = rot a (right panel). The results are given for the cylindrical symmetry with two charges located 
at the points with z =0 and 1. 
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Here the brackets denote the scalar product of charge 
vectors and 5 is the delta-function source of unit intensity. 
Then Eqs. achieve the following form : 



P = g ipi QxP 
Q ^ g ip2 P xQ 



(20) 



This system of equations (1^0)) describes the rotation of 

— >^ — + 

color charges w.r.t. to the constant vector fl —'■Pi P+ ^2 
Q with the frequency g\i~i\. The matrix C is conserved 
dC/dt = 0, so the potentials of the scalar (p and vector 
field a are quasi-static, indeed. 

The set of equations (fT9|l has a transparent physical 
interpretation. The color field - generated by two point- 
like sources - is the source of a color charge itself. Thus, a 
self-consistent environment of color charges and the cor- 
responding currents is established in the space around the 
charges. The solutions of this system were investigated 
in detail both analytically and numerically in Ref. jioj . 
There, it was demonstrated that the interaction between 
color charges is Coulomb-like if the coupling constant is 
not large, i.e. 5'^/(47r) < v^. The corresponding solu- 
tion can be well approximated by the Coulomb poten- 
tials (pi,2 = l/(47r|a; — a;i.2|) with the vector field a that 
is generated by these potentials in the next iteration in 
g. This is illustrated in Fig. [1] by drawing the "lines 
of constant force" for a (left panel) and equipotentials 
for the induced chromomagnetic field H^p = rot a (right 
panel), respectively. A line of force is a directed curve in 
the color field such that a forward tangent at any point 
shows the direction of the chromofield intensity. It is 
seen that the vector field of two color charges looks like 
the field of a permanent magnet with poles being placed 



at the points Xi,X2- The vector field a develops a con- 
stant longitudinal component |a||| — 1/{4:Tt\xi — X2\) only 
along the straight line connecting the sources at zi = 
and Z2 = 1, as demonstrated in Fig. [1] where the cur- 
rent lines of the vector field a are presented. The field 
drops sharply perpendicular to the line on which the color 
charges are located. The chromomagnetic field includes 
a single (vortex) component and, for example, for color 
charges at the points and 1 on the axis z, is given by 
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1 /I 



(47r)2 r 
•Px Q; 



1/2 



HipPxQ 
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1 



(l-2r^ + r2)i/2 

(21) 



where r, fj, are the spherical coordinates measured from 
the origin and fi = cos 9 with the angle 9 measured w.r.t. 
the z axis. The contour lines of the H^p component of the 
chromomagnetic field are presented in the right panel of 
Fig. [T] It is noteworthy that the lines of constant force 
and equipotentials seem to be mutually orthogonal at 
each z, r point. Equation ([2T|) becomes more complicated 
for other choices of the charge location. 

The total energy concentrated in the color field may 
be estimated as 



E = d 



E +H 



\Xi - X2 



Xi - X2 



(22) 



with ag = 5^/(47r), / = (6 - n^/2)/A; the terms of 
the source self-interaction are not included here. These 
estimates of the energy density of the chromomagnetic 
and chromoelectric fields clearly demonstrate the rela- 
tive contributions of the first and subsequent iterations as 
well as the dominant role of the first iteration in the cou- 
pling constant relevant for actual applications in ultra- 
relativistic nucleus-nucleus collisions. 
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A closer analysis shows that in the locations of the 
charges the potentials ipi have a singularity. The limiting 

values Vi depend on the path by which one arrives at the 
charge position. For example, for one of the charges in 
the origin we get 



1 1 
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where a is the value of the vector field a|| at the location 
point of the charge. Formally, the presence of singulari- 
ties leads to a contradiction with the initial assumptions 
of deriving Eqs. (PH)) because now a mismatch of the color 
vector rotation on the charge itself is possible. In fact, it 
signals the necessity of a more accurate fixing of bound- 
ary conditions on the charge (this problem is akin to simi- 
lar problems of diffraction theory). In particular, if these 
conditions are formulated in such a way that the vec- 
tor field a does not penetrate into the sources, one can 
avoid these difficulties [53|. However, the approximate 
solutions are still quite appropriate to be used for the 
coupling constants which are of interest for applications. 
Then it can be shown that the solutions of Eq. ([T^ with 
a boundary condition that leads to expelling the vector 
field from the charge should be qualitatively restructured 
when the coupling constant reaches the threshold values 



n 



g 



(Any 



(23) 



where / = 0, 1, 2, . . . are integers (of angular momentum) 
and / — corresponds to the situation with the penetrat- 
ing vector field. In addition to these critical values there 
is another one: It results from the same characteristic 
equation that we have used to obtain the limit (|23|) and 
is associated with the singularity power of (|19l) solutions 

in the neighborhood of the origin, where the source of 

— * 
the color charge P is located, (ipi— Vi,a) ^ , 

Anf" 



1^((2/- 
2 
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The condition to treat only the real- valued solutions 

n<(; + 1/2)2 (24) 

leads to the constraint widely discussed some time ago 
[sif . It was asserted that when reaching such a threshold 
the non-Abelian fields develop instability and, as a con- 
sequence, the color charge becomes completely screened. 
However, one can notice that the critical values are 
always reached before the limit ((24| comes effective and, 
therefore, a rearrangement of vector field takes place be- 
fore the complete color charge screening develops. The 
proper time-consuming calculations are possible only by 
applying numerical schemes with the vanishing diver- 
gence. Fortunately, such effects occur for coupling con- 
stants which apparently are of little interest for the ap- 
plications. 



The scalar component of Eqs. (1191) can be written in 
the form 

VE = 5 - aCE , 

where a column of two components of the chromoelec- 
tric field E = has been introduced. The last term 
includes the vector field a and can be interpreted as a 
charge density (column of charges Gi, G2) of the gluon 
field 



G 



-aCE 



Intuitively, it appears plausible that if the color sources 
are antiparallel, then the charge of the non-Abelian field 
could be zero. Indeed, if P = — Q, then by calculating 
the charge density components Gi = G2 = aV(<I>2 — 'I'l), 
it becomes clear that the charge of the non-Abelian field 
indeed vanishes G = G1P + G2Q = ()■ In contrast, we ex- 
pect that the non-Abelian field charge becomes maximal 
if the color charges of the particles are parallel P — Q. 
In this situation Gi = — G2 — aV(<i>i -I- $2), where, con- 
trary to expectations, the non-Abelian field charge gets 
zero G = again. Note, however, that the charge of the 
non-Abelian field will be nontrivial for other configura- 
tions. 

It is important to notice that the energy of the system 
of the non-Abelian field and the particles can be written 
(if based on the solutions of Eqs. in the following 

form 



E 



~2 ~2 

E +H 



(fx (-LP S + ip G 



(25) 



where in accordance with Eq. ([TTl) S — SiP + S2Q, i.e. 
the chromomagnetic component effectively looks like it 
would be dropped. This form reminds us of the well- 
known expression of electrodynamics for the electrostatic 
energy of a system of charges 



E 



where pi is the potential (that can already be consid- 
ered as an observed quantity) at the location point of 
the charge e^. Comparing these expressions we are able 
to conclude clearly about the role of the non-Abelian 
charge cloud and, taking into account the sign of its con- 
tribution to the energy (see Eq. (1^^ . for example) we 
find that an extra repulsion is developed in the system 
of non-Abelian charges. We speculate that this extra re- 
pulsion might also lead to an enhancement of the elliptic 
fiow V2 of hadrons in relativistic heavy-ion collisions that 
has not been considered so far. 

In view of the classical field dynamics we realize that 
Eq. (j25|) serves as the only source of information on the 
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forces acting in the system. It is interesting to remark 
that e.g. the Wong force (invented to describe the dy- 
namics of color charges in an external field) is associated 
only with the first term ~ ip5 of Eq. p5|) . as it should be 
in the chromostatic case, since it is the only combination 
permitted to construct properly its covariant momentum 
extension. It can also be constructed by counting the 
Wong force for the first and second charge explicitly (us- 
ing the definition of the chromoelectric field strength E) 

Fp = PE(xi) = V(pi\x, + V^2\x, {PQ) , 
Fq = QEix2) = V^ilx, (PQ) + Vip2\x, 0" ■ 

Singular contributions describe the self-interactions and 
should be regularized as in electrodynamics. Equa- 
tion (|25p demonstrates explicitly that the superposition 
principle for a color charge is non-applicable in the gen- 
eral situation and one has to accurately account for the 
non-Abelian field cloud energy. The structure of the 
cloud environment and, hence, the system energy de- 
pends, generally speaking, on the choice of the boundary 
condition for the charge. 

It is relevant to notice that the iterative perturbation 
series can be analyzed (when dealing with the Yang-Mills 
systems ([5])) only when discarding the compatibility con- 
dition (fTO|) . The correct construction of the iterative se- 
ries inevitably leads to the conclusion about the non- 
conservation of classical color charge. This property is 
essentially different from naively expected small correc- 
tions to the Abelian solutions because of the nonlinearity 
of the Yang-Mills equations. In a sense this conclusion 
can be deduced from the inverse statement that the solu- 
tions are automatically correct with accuracy of order g 
(in reality g"^, sec ((25|) ) if the compatibility conditions are 
taken into account. In the non-Abelian case it is tech- 
nically difficult to operate with the different gauges and 
besides, the set of physically justified (gauge invariant) 
observables is less than in electrodynamics. 

Since the strengths of chromoelectric and chromomag- 
netic fields are not proper observables, the problem of in- 
terpreting their solutions arises. Actually, in the example 
discussed above the gauge is implicitly fixed by selecting 
a special type of solutions. If one makes a choice in favor 
of the rest frame of the rotating charges, the potentials 
<Pi,2 approach some constants at spatial infinity. Such a 
gauge is inconvenient to be used in numerical calculations 
and our experiences suggest that the Coulomb solutions 
are most informative. 

Adding even one more color charge makes it impossi- 
ble to investigate the problem at the same level of consis- 
tency as for the two-particle problem because the com- 
patibility equations (jl9p cannot be integrated anymore. 
Nevertheless, it seems reasonable to consider that the sys- 
tem energy is given by the same expression (1251) with the 
Coulomb interaction of the components. The charge con- 
tributions of non-Abelian fields (diffusive color clouds) 
are small for the coupling constants of actual interest be- 



cause in the relativistic heavy-ion situation it is impos- 
sible to factorize the equations in a general form and a 
heavy problem how to construct an approximate solution 
raises. 



IV. THE FIELD OF TWO COLOR CHARGES 

Now let us consider the field of two color charges P, Q 
moving along the z-axis towards each other with veloci- 
ties V and — u", respectively. If their encounter is supposed 
to take place at time i = 0, their coordinates in the labo- 
ratory system are given as z\ = vt, zi = —wt. Obviously, 
when the particles are far from each other their inter- 
action can be assumed weak and their charges remain 
constant with high accuracy. 

Adapting the typical scale of non-Abelian (strong) in- 
teractions of about 1 fm we assume that the interaction 
is switched on just at this distance (scale) denoted as D. 
Obviously, this generates the characteristic time scale in 
the problem. Then the first "milestone" appears when 
the particles enter the interaction area, i.e. 

T = — . (26) 

The second time instant of importance is t'-^ when a signal 
of the appearance of charge Qt at the distance of 1 fm 
reaches the first particle (with the color charge P). Sim- 
ilarly for the second particle, the time when the charge 
Pt comes at the same distance is denoted by t^-, where 
Pt, Qt are the color charges before entering the zone of 
interaction, i.e. before the time T, 

. \ — W ,1 — 17 , , 

A = ^—T, t'^ = -—T. (27) 



1 



w 



These expressions are easily extracted from the scheme 
in Fig.[2]by writing out the relations for the arrival of the 
"light" signals and charges at the points of interest. In 
the figure, these events take the form of the correspond- 
ing triangles. Just at these times the color charges start 
to rotate with respect to their constant color vector that 
is the vector peculiar to the partner charge before it en- 
tered the interaction area. This regime is going on up to 
the time 

.".i-^l-^T, (28) 
1 + t; i + i« 

which is the same for both particles. It is also the time 
necessary for a signal to reach the partner providing in- 
formation on the beginning of rotation (see Eqs. ((20)) ) 
starting from its asymptotic charge value Pt, (Qt)- It 
is noteworthy that the same approach is used for formu- 
lating properly the Cauchy problem for retarded equa- 
tions [52|. 

The velocity of a relativistic particle is determined by 
the relation 

1/2 2 

1 - ^ , (29) 



TO 
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Figure 2: The time scheme for the meeting of two color par- 
ticles. The solid lines are the particle trajectories, the dashed 
lines show projections on the time axis, the dotted lines are 
the light signals. Here T is the initial time of the interaction 
when the particles become closer than the distance D = 1 fm; 
t'l is the arrival time of the signal to the first particle about 
the presence of the color charge Qt at 1 fm; similarly for the 
second particle; t'2 is the time when the charge Pt is seen at 
the same distance, where Pt, Qt are the charges before en- 
tering the interaction zone, i.e. before the time T. Until the 
time t" , the first and second particles rotate with respect to 
constant color vectors Qt, Pt- 

where m is the particle mass and £ its energy. It allows us 
to estimate the order of magnitude for the characteristic 
time in the problem as 



The collision energies of present heavy-ions facilities 
(RHIC and LHC) allow us to estimate the corresponding 
factors as f /m ~ 10 — 10^ and more. Another interesting 
scale is the distance S between the dipole charges and we 
assume it to be of the order of the interquark distance in 
the nucleon; i.e. we use 6 = 1 fm. Then in the laboratory 
system the dipole size (due to Lorentz contraction) will 
be 

Let us denote the time corresponding to this scale as 
is ~ m/£ (see below). Thus, for the relativistic prob- 
lem of interest we obtain the following time hierarchy of 
interaction stages 

t" < < ^3 < T 

(see Appendix ). 

Due to the chosen geometry of the problem an approx- 
imate solution of the Yang-Mills equations for two color 




t' 



1 2 

Figure 3: The scheme for determining the retardation time for 
the collision of two charges. Here zi and 22 are the current 
positions of the particles at two certain time moments. The 
charge of the first particle is characterized by the retarded 
time t** as well as by a segment D which measures the dis- 
tance between the particles for this time. The parameters t* 
and d for the second particle have a similar meaning (see also 
the text). 



charges can be represented as the following superposition 
(see also ((Til)) : 

^ = [^lP]t' + y2Q]t' , A,=V [^lP]t'~W [^2Q]t' ■ (30) 

The scalar potentials ^pi and (p2 can be specified as: 

^'^ 47r i?i - w i?i ' '^^ 47r i?2 - w i?2 ' 

The factor g is taken out since the normalization of the 
color charge vectors of the particles is chosen as unity. 
The consistency conditions (ITT)) include the potentials at 
the location point of the particles. Let zi and Z2 be the 
coordinates of the particles at some moment (see Fig. [3]). 
From this picture one can conclude that the potential 
values of interest are determined by the distance (it is 
convenient for a presentation in formulae to put here the 
observation point Tq at the origin of coordinates) 

Ri = d + Z2- Zi=t* , R2 = D + Z2-Zi= t** 

with d = vt* and D = wt** . Using these relations, we 
obtain 

1 9 1 9 

47r Z2 — Z\ AtT Z2 — Zl 

Thus, the Lorentz factors are compensated in such a way 
that the scalar potentials at the location of the charges 
can be simply expressed in terms of the distance between 
the particles in the current moment. Also, taking into 
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account the form of the vector potential, we arrive at the 
following expression for the compatibility conditions : 



P = a, Q{t-tll) X P , 



Q = a 



' \Z1-Z2\ 
\ + V w 
\Z1 - Z2\ 



(32) 



The retarded time can be obtained (using Fig. [3l for 
example), as 



t% 



w 



1 



t 



,** 

'•12 



1 — W 



t , for i < 



For t > (in the above formulae) the following replace- 
ment should be made: v,w ^ —v,—w. In the electro- 
magnetic case, the charge is conserved and its retarded 
time is defined simply as t — t' + R{t')/c [s^. Below 
we will apply a reduced notation for the retarded time 
without any further explanation. 

It is convenient to introduce the variable = — In 
by means of which Eqs. p2l) are reduced to 



(33) 



F' =w Q(V3- A*) X F , (t<0) 
Q' = P((p - A**) X Q , 

for negative times, and 

P' = ~Lo Q{ip + A**) X P , {t>Q) 
Q' = -uj F(^ + A*)xQ, 



(34) 



for positive times. Here the prime denotes differentiation 
with respect to Lp, 



1 



V + w 



A* =ln 



1 



1 



A** = In 



1 



w 



1 



The procedure for obtaining an approximate solution 
is given in the Appendix. To find the strength of the 
chromoelectric and chromomagnetic fields, it is also re- 
quired to know the derivatives of the color charge vectors 
w.r.t. the retarded time dP{t')/dt', dQ{t')/dt'. They are 
computed by the explicit formulae ([32t . 

As an example, let us consider the field created by 
two relativistic particles moving with velocities: v ~ 
1 — 2 -10"^ and \w\ = 1 — 1-10~^. The energy-mass 
ratio for the first particle is S/m ~ 16, and for the 
second particle is 8/m. ~ 22. Let us take the coor- 
dinates of the observation point as rp — (2,0,1), i.e. 
X = 2 fm, z — 1 fm. For comparison, we consider 
also the field created by particles with electric charge 
±e of the same interaction strength as the color charges 
(e^/(47r) = /{A-k) — 0.3) moving with the same veloc- 
ities. The potentials for the scalar Lp and z component 
of the vector fields of color particles are presented 
in Fig. 21 where the initial angles in the color space are 
determined as 6* = 7r/1.95, </) = 7r/20 for the first par- 
ticle and 6 — — 7r/1.95, 4) = — 7r/20 for the second one. 



In this case the color charge of the first particle, for ex- 
ample, will be P = (Pi,P2,P3), with Pi = sin cos i^, 
P2 = sin 6* sin 0, P3 ~ cos 9. This configuration of color 
charges P, Q corresponds to almost oppositely directed 
charges at the initial stage. The dashed lines in this figure 
show the three color components of the scalar potential 
(left panel)) and the third component of the vector po- 
tential (right panel) . The open circles are plotted for the 
potentials corresponding to the enhanced electrodynamic 
coupling. The modulus of the scalar \ip\ and vector \Az\ 
potential are displayed by the solid lines in Fig. 01 We 
recall again that these variables are not directly observ- 
able. The selected observation point rp = (2, 0, 1) intro- 
duces some asymmetry reminding of the asymmetry in 
the transverse plane for peripheral nucleus- nucleus colli- 
sions which lead to some dominant component (33 |. Ac- 
cording to the choice of the geometry, we see that the first 
maximum corresponds to the passage of the first particle 
at the closest distance to the observation point and then 
the second maximum is the corresponding passage of the 
second charge. Both bumps are located symmetrically 
with respect to the point t — 0. For the selected con- 
figuration of color charges, one of the components of the 
potentials - shown by the dashed lines - dominates and 
almost coincides with the appropriate modulus of vectors 
in the color space (solid line). The meeting point is lo- 
cated at a distance of P = (x^ -|-z^)^/^ = \^ ~ 2.24 from 
the observation point. It is seen that there is a notice- 
able difference between the scalar and vector potentials 
as compared to the electrodynamic case (open circles) at 
the appropriate time when a signal on a meeting of the 
particles has arrived. As discussed in the Appendix, the 
rotation of color charges is described by Eqs. (p9)) -(|43 l) . 
The effective rotation frequency in the neighborhood of 
the meeting point t is lu" = aq/ln\t\, i.e. the color 
charges are rotating infinitely fast. 

The time dependence of the strength components of 
the chromoelectric (left panel) and chromomagnetic 
Hy (right panel) fields is shown in Fig. [5] by the three 
dashed hnes; the sohd lines correspond to the modulus 
of the chromofields | \ and | Hy \ . We have cut the sin- 
gular peaks at some threshold and thus the lines look 
somewhat irregular. In full agreement with Eqs. (jlSp . 
the Ex and Hy components are dominating. In both 
cases two maxima (minima) caused by passing the color 
charges in the vicinity of the observation point are clearly 
visible. Note that here (and in what follows) the color 
field strength is plotted in dimensionless units where c 
is the color charge of the appropriate field component 
and m^r is the pion mass. The charge velocities con- 
sidered roughly correspond to the RHIC energy where 
the maximal electromagnetic field eHy/rn^ reaches a few 
units (3JI . This value is essentially smaller than those in 



the color charge case (see Fig. [5]). It is seen that for the 
color charge configuration considered the chromoelectric 
and chromomagnetic field are quite similar to the field in 
the case of enhanced electrodynamics. Some difference in 
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-10 12 -10 12 

t Ifm/c] ^ t [fm/c] 

Figure 4: (Color online) Time evolution of three color components for the scalar <j) (left panel) and z component of the vector 
potential Az (right panel) for two oppositely directed charges in color space as a function of time is given by the dashed 
lines. The solid line shows the absolute value of the isovector potential \<j)\. The electrodynamic potential with charges ±e 
corresponding to the coupling constant Ue = 0.3 is displayed by open circles. 
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Figure 5: (Color online) Three color components of the chromoelectric E (left panel) and chromomagnetic H field strength 
(right panel) for two moving color charges of the opposite signs as a function of time are given by the dotted lines. The solid 
line shows the absolute value of the isovector potential [E[. The open circles show the electrodynamic vector field corresponding 
to charges ±e and coupling constant Qe = 0.3. 



the hight of the first two maxima are caused by different 
velocities of color charges. A significant difference at the 
third maximum is due to the arrival of a signal from the 
meeting point of particles to the observation point, where 
there is a noticeable additional contribution of chromo- 
electric and chromomagnetic fields associated with the 
temporal change of the particle color charges. 

The modulus of transverse = \/^ ^'^^ longitu- 
dinal I-Ezi = component of the chromoelectric field 
is shown in Fig. [6] for two configurations of the initial 
color charges. The longitudinal component is strongly 
suppressed, as it should be due to relativistic effects, 
but the signal from the meeting point of the particles 



leads to almost equal contributions. It is noteworthy that 
the change of the initial configuration of color charges 
from parallel to orthogonal (cf. left and right panels in 
Fig. [6]) does not change the evolution in the absolute val- 
ues of the chromoelectric field though the sign of the color 
Ex,Ez components are different and change with time. 
It is also important to note that every collision noticeably 
changes the position of color charges in the color space. 

Characteristics for the case of initially almost orthogo- 
nal color charge vectors are presented in detail in Figs. [7] 
and [51 These results are obtained for particles with color 
charges defined by the following angles in color space: 
6» = 0, = 7r/20 for the first particle 9 = -7r/1.95, 
(f) = — 7r/20 for the second one. We shall henceforth refer 
to this case as non-Abelian. The notation in these figures 
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Figure 7: (Color online) Non-Abelian configuration of color charges (see text). The three components of <j!> (left panel) and Az 
(right panel) are shown for the initially orthogonal charge vectors. The notation is the same as in Fig. |4] 




-10 12 -10 12 

t [frn/c] t [/m/c] 



Figure 8: (Color online) Non-Abelian configuration of color charges (see text). The three components of chromoelectric E (left 
panel) and chromomagnetic H (right panel) fields are shown for the initially orthogonal charge vectors. The notation is the 
same as in Fig. [S] 
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is identical to that in Figs. |4] and [5] The main result here 
is that in this case there is no dominant component, as in 
the case of mutually opposite charges considered above, 
but two preferred directions in color space are significant, 
as evidenced by the corresponding maxima shown by the 
dashed lines in Figs. [7] and El 



V. THE FIELD OF A COLOR CHARGE AND A 
COLOR DIPOLE 

As in the previous Section, let us assume that there is 
a particle with color charge P which moves along the z 
axis with velocity v. We denote it as the first particle. 
Suppose that the dipole made up by the second and third 
particle with charges Q2 and is aligned along the z 
axis and moves in the opposite direction with velocity w. 
We denote the distance between the charges in the dipole 
as 5w in its rest frame which is taken as 1 fm. In the lab- 
oratory frame the charges are located closer to each other 
= (1 — 'uP'Y^''^5yj due to the Lorentz contraction. Let 
the first and second particle encounter at zero time be- 
5' 

fore the time t-^ — ^ of the second meeting between 
the first and third particle (see Fig. [9]). Then the trajec- 
tories of the first, second and third charges are defined as 
Zi = vt, Z2 = —wt, Z3 = —wt + 5!^. Similarly to Eq. (1501) . 
let us adopt an approximate solution of the Yang-Mills 
equations in the form 



If = [(piP]t' + [f2Q2]t' + [fsQslf , 



(35) 



The potentials at the location point of charges can be 
found in the same way as in the case of two color charges; 
this leads to the compatibility conditions for each of the 
charges, see Eq. (ITU) 



1 + V w 



^ \Z1 - Z2\ 

Q3{t~tXT)-KP , 



Q2{t-tll) xP + ag 
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a. 
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' \ZI-Z2\ 

Qsit - t*2T) X Q2 



P{t-t;^) xQ2 + ag 
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(36) 
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\zi - zal 

Q2{t-t**2)xQ3 , 

With tl, = -\^{t - t,), = -j^{t - h) for 
t3 < t and q* = ts, til* = for t < h 

(with the substitution v^w — >■ ~v,—w). The retardation 
times are determined similarly to the case of two color 
charges, where, in particular, the retarded times f^i and 
are given (see the relevant scheme in Fig. [3]). An 
interesting peculiarity of the resulting system of equa- 
tions (1361) is a strong suppression of the contributions of 




Figure 9: Scheme of the meeting of a color particle and dipole. 
The particle trajectories are plotted by the solid lines, the 
dotted lines are light-cone lines corresponding to the arrival 
of the signal from the particle meeting point. The arrival of 
the signal is marked by the points 12 and 13 (see text). The 
time is the moment when the first and the third particle 
meet each other. 



charges 2 and 3 flying in the same direction, which enter 
into the equation with the Lorentz factor 1 — . Thus, 
in the ultra-relativistic case the mutual influence of color 
charges flying in the same direction may be ignored and 
the system can be considered as "frozen" . This is true 
even in the cases where the system has the size of a nu- 
cleus. 

We are interested now in the particular case when the 
color dipole charges in the initial state are opposite, i.e. 

— ~Q2- The prescription for obtaining an approxi- 
mate solution in this case is given in the Appendix. 

The strength of the chromoelectric and chromomag- 
netic fields - generated by particles moving with the ve- 
locity V — 1 — 2 10~^ and by the dipole with the velocity 
|w| = 1 — 1 10~^ - is demonstrated in Fig. [TOl The color 
charges of the particles at the initial time are determined 
by the following angles in the color space: 9 = 7r/1.95, 
(f) = 7r/20 for the first particle, 6* = 0, = — 7r/20 for the 
second one and with the opposite angles for the third par- 
ticle. This configuration is denoted as 'non-Abelian'. We 
do not show explicit data for the potentials, as they pro- 
vide little information. It is seen that generally two first 
maxima of the field strength generated by color charges 
in the vicinity of the observation point are reasonably 
reproduced by the Coulomb-likc solution (dotted lines in 
Fig, nop . However, the first maximum formed in passing 
of the dipole is not smooth but has an "up-down jump" 
(or "zig-zag") shape. This structure is caused by the 
opposite color charges forming a dipole. 

The observation point is located now in another place 
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X = 1, z = 1 fm as compared to the case of a meeting 
of two particles at a; = 2, so the signal of the meeting of 
the color charge and the dipole arrives at a time about 
r = (x^ + z^)^/^/c = \f2 ~ 1.41. Here, at the meeting 
point (r = 1.41) the Coulomb-like solution predicts no- 
ticeably narrower distributions due to the so-called " color 
charge glow" which results from the color charge interac- 
tion through the time dependence of the color vectors D 
(the last term in Eqs. ([T5|) ). 



VI. THE FIELD OF TWO COLOR DIPOLES 

Let us finally consider the field created by two color 
dipoles. The first dipole is formed by the first and fourth 
particle with charges P\ and P4, respectively. The second 
dipole is made up of the second and third particle with 
charges Qi and Qa. At the initial time the color charge of 
each dipole is neutral: Pi = — P4, Q2 = —Qa- Particles 
of the first aligned dipole move along the z axis with the 
velocity v and the second aligned dipole moves towards 
them with the velocity w. Let us denote the distance be- 
tween the charges in the first and second dipole in their 
rest system as Sy and S^, respectively. In the laboratory 
frame, these distances are contracted (5(, = (1 — v-^Y^^Sv 
and S'^ — {1 — w'^y^^Sw Again, it is convenient to intro- 
duce the time scales 

, _ S'^ K 

13 — — \ , 14 — — ; , 

V + W V + w 

when the first particle meets the third one and the fourth 
meets the second one, respectively (see Fig. [TT|). At the 
time ^3 -I- t4 the fourth particle meets the third one. 
Suppose for convenience that the meeting of the first and 
second particle occurs at time zero. Then the trajectories 
of charges with numbers one, two, three and four are 
defined as zi = vt, Z2 = —wt, Z3 = —wt + 5'^, Z4 = 




4 1 ? 3 



Figure 11: Scheme of the meetings of two dipoles. The solid 
lines are the particle trajectories. The meetings points of the 
particles are marked by ta for the particles 1-3, t4 for the 
particles 4-2 and ts + for the particles 4-3. 



vt — 5'y. In this case the ansatz for the superposition of 
approximate solutions has the form 

if = [(piPi]t' + [mQ2]t' + [vaQs]*' + ^4^4]*' , 

.(37) 

^ V [(flPl]t' - W [ip2Q2]t' - W [(ySsQs]*' + V i^iPilt' ■ 

Similarly to the previous Sections the compatibility con- 
dition (see Eq. for each charge can be written as 
follows : 

~ 1 + vw ~ , ~ 1 + vw 

Pi = ag r Q2(t - ti2) X Pi + ag r 

W - Z2\ \zx - Z3I 
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1^1 — ^3 



?c /, 5 . \ ~ ~ jy ~ substitution should be made — > — u, — w); tl?* — 

Q,it - ) X + P4(t - t,, ) X P„ -^(t - t3 - ^,), = - t3 - ^4) for 

~_ 1 + vw ~ ~ 1 — w'^ i < t3+t4 (for i3+i4 < t, the substitution should be made 

O X O + a ^ ^ ^ ^ P (t t^^) X O ^^^'^ Appendix we present the procedure for obtaining an 

23 ) V2 ff — Z4I ^ approximate solution for two color dipoles. 

138) 

1 + w w ~ ^ ~ 1 ~ w In Fig. [T21 the strength of chromoelectric and chromo- 

^3 + 0!g _ magnetic fields created by two dipoles with the veloci- 

1 -(_ ^ _ _ ties V = 1 — 2 10^^ and |w| = 1 — 1 10^^ is demon- 

■ ^2(^ — ^32) X (33 + Qfg 1 r P4{t ~ t^^) X Q3, strated. The color charges of the particles at the ini- 

'^^ tial time are determined by the following angles in the 

^ / " P^(t _ t*) xP^ + a, 1±1^ color space: 9 = ^/1.95, = n/20 for the first par- 

|zi — Z4I 1^4 — 2:21 tide and 9 — 0, (p — — 7r/20 for the second one and 

~ ~ I ^ y ^ ^ with the opposite angles for the third and fourth parti- 

■ Q2(i - tl^) X P4 + Qfg _ Qsit - til*) X Pa, cle. As is seen (by the dotted lines in Fig. the two 

first maxima are reproduced by the Coulomb-like solu- 

, ,.rr , in 1 ,** v + Ti!/, tlou includiug thc up-dowu lump cffcct for both maxlma, 

where the different times are defined as t49 = ^ (t— , , , • r. -,r r^i 1 ,• ■ . ^ ,i 

1 — w ^ noted above m Section V. ihe observation point has the 

^4)5 tiX — ~ 1 ^ ^ ~ ^4) for t < t4 { for t4 < t, the coordinates x = 1 fm and z = 1 fm, so the signal on the 
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meeting of the color dipoles arrives at the time about 
T = (x^ + z2)i/2/c = 72 ^ 1.41. In fact, the whole 
area from the second maximum till the meeting point is 
filled by the chromoelectric and chromomagnetic fields of 
a noticeable strength and near t — 1.41 this looks like a 
band. To clarify the band structure, we zoom the region 
near the meeting point in Fig. 1131 It is seen that really 
this area consists of four maxima of huge intensity which 
in dynamical systems can result in large local color fluc- 
tuations. This shining of the color charge glow is mainly 
due to color charge interactions through the time depen- 
dence of the color vectors D (the last term in Eqs. 
), which is certainly not reproduced by the Coulomb so- 
lution with the enhanced coupling. It is noteworthy that 
the color charge glow effect is getting even larger for more 
complicated systems. 



VII. CONCLUSIONS 

In this study we have considered several elementary 
configurations of relativistic partons (which are charac- 
teristic for the initial state of ultra-relativistic nuclear col- 
lisions at the quark level) with the non-Abelian charges 
for the SU{2) group in the classical limit. It is shown 
that as in the case of non-relativistic particles, the system 
generally shows Coulomb-like features, and the analogy is 
convincingly supported by a comparison with electrody- 
namics for the same coupling strength. In distinction, we 
find an additional strength of chromoelectric and chro- 
momagnetic fields close to the meeting of particles which 
is caused by the explicit time dependence of the color 
charge vectors. In the chosen gauge the interaction of 
the color charges results in the rotation of the color vec- 
tors which becomes very fast close to the meeting point of 
two particles. The chromo fields are essentially stronger 
than the corresponding electromagnetic fields created by 
a moving source at similar conditions. The longitudinal 
and transverse field components of this signal are of the 
same order of magnitude in contrast to the longitudi- 
nally compressed field (due to the Lorentz contraction). 
Changing the observation point one can see the predicted 
shift of the meeting signal of the particles. Really it turns 
out that for complex systems the whole observation zone 
of a reasonable size of an order of a few fm near the meet- 
ing point is filled up by this intensive signal of a strong 
interaction scale which we have denoted as " color charge 
glow" . 

In fact, there is no direct evidence of the Coulomb law 
for the interaction between quarks. The leading hypoth- 
esis to explain the observed behavior of quarks is the 
idea of vacuum gluon fields. In this paper we have not 
touched on this topic but one should note that estimates 
based on the instanton vacuum model indicate that the 
vacuum field may be much stronger than the fields gen- 
erated by the collision of quarks ;53| . Every open color in 
the instanton medium is screened which indicates the im- 
possibility of forming gluon fields of noticeable intensity 



and, in physical meaning, apparently implies the trans- 
formation of the gluon field quanta in energetically more 
favorable configurations. 

Using the proposed approximation that particles be- 
gin to feel the presence of the non-Abelian charge of a 
collision partner only below the distance of 1 fm, the dis- 
cussed color configurations can apparently be considered 
as a model for a color pre-quark-gluon matter. As to a 
possible formation of some objects like a color-flux tube, 
a larger system similar to that formed in ultra-relativistic 
heavy-ion collisions should be considered. Here new col- 
lective effects like the Debye screening come into play, 
additionally. However, these effects are beyond the aim 
of our study and are the subject of future investigations. 
Also, an important "practical" conclusion is that, in 
essence, the correct assessment for the strengths of chro- 
moelectric and chromomagnetic fields is obtained just in 
the approximation of pure "electrodynamics" with en- 
hanced coupling (7". At the same time, the experience 
gained so far [5J| indicates that such a signal of "color 
charge glow" might manifest itself not so directly as sug- 
gested by the figures above. Although only basic con- 
figurations of moving partons have been examined, the 
qualitative picture obtained should hold also for more 
complex color configurations. We speculate that signals 
of the additional repulsive non-Abelian interaction might 
leave its traces in the early development of the elliptic 
fiow in relativistic nucleus-nucleus collisions prior to par- 
tonic equilibration. 
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Appendix 

Time hierarchy of interaction stages 

Below, in the discussion of a color charge, a color dipole 
and two dipoles we need to compare the description at 
different scales. It is convenient to specify them by the 
following items: 

1. The time t <T,T <t <t" (recaU t" < 0) 

The form of the solution for negative times has al- 
ready been specified above, for example, the color 
charge vectors are constant up to the scale t' 

P = Pt , t<t[ 

^ ^ (39) 

Q = Qt : t<t'^. 

Further, up to the time t" , each charge rotates rel- 
ative to the constant vector of the opposite charge 
particle. Let us introduce the basis (static) vectors 
in the color (isotopic) space, through which it will 
be convenient to express the solutions of equations. 
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The triple orthogonal unit vectors are associated 
with the vector P in the form: 



Qt 



Pt X 



where cos 6' = [PtQt)- A similar basis associated 
with the vector Q is defined as 



-np„ 



Pt X nq^ 



The solution of Eqs. for the times considered 
can be represented as follows 



P< = cosO Qt + sin 6* {cos[uj{ip — ip'i)] mp^ 
sm[uj{ip — ip[)] np^} , t[ <t < t" , 

(5< = COS0 Pt + sin 6* {cos[uj{ip — (/Jj)] fiiQ^ 
s\n[uj{ip ~ ip'^)] UQ^} , t'^<t <t" , 



(40) 



where we use the notation Lp'-^ = —hi\t'i\ 
-\n\t'\. 



V'2 



The time t" <t< -t" 

Now consider the asymptotic solutions for large (f 
(in the neighborhood of the meeting point of par- 
ticles). In this case, one can neglect the retarda- 
tion factors A*, A** in Eqs. ([331), ([Ml)- Then, for 
both positive and negative times these equations 
describe the rotation of the color charge vector with 
respect to the vector P + Q which is conserved. It 
seems reasonable to approximate the solutions of 
Eqs. ([5^ on the whole semi- axis of the negative 
time including the time of the meeting of particles, 
by matching the solutions on the t' scale with those 
on the t" scale, where the retardation is neglected. 
The numerical study of the system ([5^ justifies in 
general such an approximation. Thus, at the time 
t" the color vectors are 

P" = P<(t") , Q" = Q<(0 . 

Therefore, the vector P-|- Q is defined with respect 
to which the vector of the color charge particles is 
rotated on the scale T" < t < 0. Similarly to the 
considered case of solutions for large time t', let us 
introduce three orthogonal unit vectors 



P" + Q" ^ P" xfl 
, npii^^ 



\P" + Q" 
rhpn = X npn , 



IP" X n\ 



by means of which the solutions of systems 
be expressed as follows 

P< = cVl + s {cos[w"((^ — ip")] fhpii 

- sin[w"(<^ - p")] npn} , 

Q< = c VL — s {cos[uj"{ip — ip")] rhpn 

- s\n[uj"{ip - ip")] npn} , 



(41) 



uj" = uj \P" + Q"l c = cos 



(P"Q") 



. (P"Q") 

s = sni 



2 

p" - -ln|t"| 



This solution describes an infinitely fast rotation of 
the charges in the neighborhood of their meeting 
point, if they are treated as a function of time, and 
we have a rotation with a constant frequency u" in 
the logarithmic variables p. To continue the solu- 
tion beyond the singularity for positive times, we 
assume that in the time interval — ^min < ^ < ^miiu 
the particle charge vector does not change. The 
time scale ij^in of artificial character and can be 
chosen arbitrarily small, in principle. Without re- 
tardation the solution of Eq. (p4l) at < > t^^^i can 
be written as 



P> = c + s {cos[a;"(iy9 — P*)] mpn 

+ sin[ijj" {p — p*)] npn} , 

Q> = c fl — s {cos[a;"(iy9 — p**)] rhpn 

+ sm[uj" {p — p**)] npn} , 



(42) 



where p — ~lnt, p* , p** are arbitrary phases to be 
defined by the known values of the charge vectors 
on the left end of the segment — imin 



-P<(~*min) 
^min) 



P>{t 



min) 



Q<{- ''-mm ) = Q>{t mm) 



Comparing Eq. and (|^ . one can obtain re- 
lations between the phase p" and phases p* and 
p** . It is seen that it is convenient to put the time 
^min such a way that sm[uj" {p-^^^ — p")] = 0, 
with i^jnin = - Intmin; i-e. <^min = 27r n/uj" + p" , 
where n is an integer number. Then one can put 
p* = p** = p" , As a result, the expression for pos- 
itive times will not be overloaded by formal phase 
shifts and at the same time an arbitrarily small time 
scale ijnin can be chosen. Thus, the solution for the 
time scale ijnin < t < \t"\ has been obtained. As a 
result, we can see that 



p<(n = p>(|t"i) = p" 

Q<{n = Q>{\t"\)^Q" 



i.e. in this approximation the phase shift is not 
observed, and the color charge vector at the exit 
from the scale t" coincides with that at its entrance. 

3. The time: -t" <t< tout; ^out < t 

We continue the approximate solution to larger 
times in such a way that a change of the color 
rotation regime on the scale tout W\ should oc- 
cur symmetrically with respect to the negative time 
and the rotation should stop when the scale t ~ \T\ 
is reached. This solution for the positive time can 
be written in the form (it is worth to mention a 
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change of the sign in this system in comparison 
with Eq. dMl)) 



P> = cos6'o Qo + sin 6*0 {cos[uj{(p — ip*)] mp^ 
+ sin[c^(^ - ^:)] npj , \t"\<t<\t'\ , 

Q> = 008 6*0 Po + sin 6*0 {cos[uj{ip ~ ipl*)] rhq^ 
- sin[a;((^ - Lp**)] ngj , 

with yet unknown basis vectors 

~ Po ^ Qo ~ ~ 

Qo , np^ = — — — , mp^ ^ Qo^ np^ , 



(43) 



X Tic 



on which the solution is spanned and with corre- 
sponding phases ipl* . Here 0086*0 = {PoQo)- 
This information should be restored using the avail- 
able P>(|t"|) Q>{\t"\) = Q". The vector 
product P" X Q" gives another condition in ad- 
dition to the relations (|43l) . Applying the vector 
algebra rules, these relations can be presented in 
matrix form 



an Po -t- ai2 Qo + ai3 "Po f"" , 

021 Po + a22 Qo + 023 "Po = Q" , 

031 Po + Qo + 033 "Po = P" X Q" 

with the coefficients 



(44) 



an = ci , ai2 = c(l - Ci) , 013 = Si , 

a2i = c(l - C2) , 022 = C2 , a23 = -S2 , 

031 = (a - /3c)/s , 032 = (/3 - Q!c)/s , 033=75 

where the following notation is used 



a 


= -C2S1 


- c(l - Ci) S2 






= C1S2 + 


c(l - C2) Sl, 




7 


= C1C2 - 


c^(l-ci)(l- 


C2) , 


c 


— COS 60 - 


= (PoQo), 




s 


— sin do , 


C\ — COS 0* , 




Si 


= sin 6** , 




'/':) , 


C2 


= cose** 


, S2 = sin6l** 





Now, if all the phases are known, then one can find 
Po, Qo using the inverse matrix A^^ . 

In order to determine the phases, we use an impor- 
tant relation for the scalar product of vectors P" 
and Q" 



(1 - Ci)(l - C2) <? + (Ci + C2 - C1C2) C - S1S2 



(P"Q") 



(45) 



A numerical analysis with the coefficients ci, ci (si, S2) 
shows that the resulting cubic equations for the cosine of 



the angle between the vectors Po and Qo has one real 
and two complex conjugate roots. The real root does not 
always satisfy the restriction |c| < 1, i.e. in general all of 
the coefficients ci, C2, (si, S2), and c should be consistent. 
The boundaries of the acceptable region are defined by 
setting c = ± 1 , then 

±(1TS152) = (P"Q") , 

In particular, for the parallel and anti-parallel vectors 
P" and Q" we have si = S2 = 0. For an approximate 
solution we can restrict ourselves by the particular case 
Sl = S2 = just as we did in choosing the phase p^^^^. 
Then for the phase one can get 

a;((^" - (/s*) = + 27r n , ujiip" - <y9*) = tt + 27r n , 

where n is an integer. Similar relations hold for the phase 
(/?**. The time scale, when the rotation around the con- 
stant charge vector of the particle-partner stops, is de- 
fined by the condition t^xA = 6""^° . The analyzed cases 
show that for the considered particle velocities (with rea- 
sonable accuracy) will fall at the scale of if the phase 
is chosen as ^j* = V*o = (/f" — tt/w. This regime is ap- 
proximately applied when mjE < 10^^'^, besides with an 
increase in energy reaching of the required scale gets bet- 
ter. In principle, nothing prevents to install the time tout 
also for moderate relativistic velocities, because at this 
scale ~ t' the change in color charges is insignificant, as 
compared to the scale t" . With this choice of the phase, 
Eq. ([45]) takes the form 



4 - 3 c = (P"Q") . 

This equation has one real root which always satisfies the 
condition |c| < 1, and two conjugated imaginary roots. 
Now, with the known coefficients ci, C2, (si, S2) and c, 
one can find the vectors of the particles that define the 
basic three vectors which are spanned over the solutions 
for positive time scales tout and |r|. The equation set 
(|^^ for particular cases interested takes the form 

-Po + 2c Qo = P" , 
2c Po-Qo = Q" ■ 
^From here we get the solutions 



Po 



1 pn 2c 



4c^ - 1 



4c^ - 1 



4c^ - 1 



4c^ - 1 



The explicit expressions provide an approximate solution 
of Eq. (1321) continuous in time which can be applied to 
the entire time axis. 



Color charge and dipole 

As before we assume that particles begin to feel the 
presence of the third-particle color charge when ap- 
proaching the distance D estimated as 1 fm. The signal 
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on the presence of the charge P arrives at the second 
particle at time iji and Ett the third one at t'^ = + ^3- 
We ignore the time difference between t'2 and t'^ (due to 
the Lorcntz contraction), i.e. the charges of the second 
and third particles do not change up to the time • With 
the same degree of accuracy, the signals from the second 
and third charges will come to the first particle at the 
time t[ . Until this point the first particle charge P is not 
changed in time. Then the charges rotate with respect 
to the constant charge vector which the particles had at 
the entrance in the interaction zone. For these times the 
compatibility conditions (I43|) become 



P^U} 



1 



1 



.1*1 \t-t3\, 

^2^^ —Pt X O2 - ^\\ —Qt X Q2 



\t\' 



t3 



(46) 



with cjii 



7-^— tPt X Os - cjii ^Qt X Q3 , 
|r — 13! 13 



„, I . From this system one can con- 

elude that due to the factor of the charges of Q2 and 

Qz can be considered as mutually opposite well away 
from the meeting point (as well as at the entrance to 
the interaction zone), in particular, until the scale i', the 
dynamics of color charges is described by a simplified 
system of the two equations 



P^LO 



1 



1 



\t\ It-hl 



QtxP , 



(47) 



Q2 — Q, Q3 — —Q- Comparing this set of equations with 
Eqs. ([5^ we get a solution in the form 



P< = cos 9 Qt + sm 9 {cos[uj [rj — rj'i)] mp^ 

- sin[a;(77 - 77^)] np^) , t'-^<t< t" , 
Q^ — cos 9 Pt + sin 9 {cos[uj{(p — ip'2)] friQ^ 

- sm[uj{if - riQ^} , t'2<t<t", 



where 



Inlil 



■ \n\t-t3\ , 



with the initial data on the scale t' . In this case, the 
dipole charges act on the color charge of the first particle 
weaker than a single color charge because the contribu- 
tions compensate each other. This scheme is equivalent 
to item 1 of the previous Section, and describes the be- 
havior of charges for times t < T, T < t < t" . 

As in the case of the two color charges let us take as 
an acceptable approximation to the exact solution the 
matching of the solution on the scale t", which neglects 
the retardation, and the above solutions for the time scale 
larger than t' . In the meeting area the charges of the 



first and second particles are described by the familiar 
equations 



P = .^02XP 



(48) 



In contrast to the singular behavior dictated by this sys- 
tem, the charge of the third particle obeys 

Q3 ~ w Pt X Q3 . 

'3 

In many applications the time scale t" is so small that 
in this interval the color vector can even be taken as 
a constant. As was mentioned in the section devoted to 
the two color charges, after the singularity point for a 
positive time t = \t"\ the charges take the same position 
in the color space as at the entrance to the singularity 
zone. Therefore, the continuation to larger positive times 
used for the two color charges holds valid in the color- 
dipole case, the color charge dynamics being described by 
the simplified system of equations (|T7|) . The initial data 
are defined by the obvious matching of the condition with 
a singular solution at the time t = \t"\, which we do not 
present here. It is clear that the same considerations 
of the behavior of color charges can be applied to the 
meeting point of the first charge with the third particle 
on the scale of t" in the neighborhood of . 



1 



- UJ 



\t-t3\ 
1 

\t-h 



O3 X P 
rPx O3 



(49) 



In its turn, the charge vector of the second particle Q2 
can be considered in this segment as constant. As in item 
2 of the previous section for the times t" < t < —t", here 
for the times t" < t — < —t", one should construct 
appropriate three basic vectors based on the continuous 
vector P{t3+t")+Q3{t3+t") and introduce an additional 
time scale tjjiin which is defined by the relation (^jj-^iji = 
2^n/wt3+t-.+(p", nt,+t" = \Pih + t") + Q3it3 + t")\. In 
this way we obtain solutions analogous to those in item 
2 but for the case of two color charges at the meeting 
points of the particles 1-2 and 1-3. 

The continuation of the solution to larger times should 
be carried out by analogy with item 3 of the previous 
Section, i.e. first, to construct the description of charges 
on the scale —t', then on the scale of t^ + t' and so on, up 
to the scale t^+t" where there is a meeting of the first and 
third particle and the behavior of charges is singular. But 
such a meticulous description apparently is not needed, if 
approximate solutions are considered. So just accept that 
the passage to the regime of rotation around constant 
vectors of color charges of particles-partners at the exit 
from the interaction zone occurs somewhere on the scale 
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ts- Such a solution of (l47l) is as follows 



Q> = 



cos 6*0 Qo + sin 60 {cos[aj(77 - r/*)] mp^ 
sin[w(r/^- «pj , \t"\ < t < Vt ' 
cos 6*0 Po + sin 6*0 {cos[a;(v3 - ipl*)] rhq^ 
sin[a;((^ - Lp**)] ngj , 



where as in the previous Section, one should determine 
the basic three vectors and phases using available infor- 
mation on the scales t", i.e. P", Q" . The corresponding 
systems of equations have the form (|44|l , (|45l) , where the 
variable substitution — )• 77 should be made for variables 
related with the vector of the first particle charge P. As 
was noted in the previous Section, the easiest version of 
the inverse problem of reconstruction of the basis vectors 
and phases would be appropriate if this choice of vari- 
ables gives si = S2 = 0. An analysis shows that for the 
problem of the particle and the dipole one can take 

With this choice of phases, going on to large positive 
times occurs somewhere on the scale of for the first 
particle and on the scale t'2 for the second one 



tol 



to2 



Xo2 = V 



27r 

TT 



As mentioned above, there is no sense to complicate the 
task by better matching the transition regime for asymp- 
totically large times. Now Eq. (|^ becomes (ci — 1, 

C2 = -1) 

c - (P"Q") . 

The system of equations (HH) in which we are particularly 
interested takes the form 

Po = P\ 

2c Po-Qo = Q" ■ 

In this way we complete the description referred to the 
item 3 for the problem of two colors charges. The whole 
procedure to obtain approximate solutions for the sys- 
tem of compatibility equations (j36p is reduced to the de- 
scription of the behavior of two basic charges, since the 
partner charge in a dipole pair can be considered as ad- 
justed. The passage for a short time on the scale t" to a 
singular rotation regime is not accompanied by a change 
in phase. 

We have missed some interesting effect of the arrival 
of the signal to the particle-partner in the dipole from 
the meeting point of two other particles. In Fig. [9l the 
marked point 12 corresponds to the light signal coming 
from the meeting of the first and second charges to the 
third color charge. In these times on the scale t" the 
charge of the third particle is described by the equation 



P X Q3 

13 



(50) 



where the charge of the first particle P is rapidly chang- 
ing in a singular manner. Unfortunately, it is difficult 
to derive analytical expressions describing the behavior 
of the color charge of the third particle, and in this pa- 
per we simply ignore this important, but short episode. 
It is also important to note in turn that a signal about 
events happened with the third charge prior to the colli- 
sion of the first charge with the third particle reaches in 
time the first particle. Then, an signal will come to the 
third particle even before the meeting of the first and the 
third particles, and so on. There is some danger that we 
are not able to control the behavior of charges at the sec- 
ond meeting point because one should trace the ladder of 
events up to the meeting point (we have not depicted in 
Fig. [S]an appropriate sequence of signals similar to those 
shown in Fig. [5]). However, it is noteworthy that these 
processes should not change too much the color charge 
phase since the first signal comes on the scale ^ m/£, 
and the second one occurs on the scale ■ t" ~ m^/E^, 
while the merging of solutions takes place on the scale 
t". 



Two color dipoles 

We discard a detailed prescription for obtaining ap- 
proximate solutions, as was done in the previous Sec- 
tions. Let us write down directly the simplified system of 
equations which allows to describe the behavior of color 
charges for large times. 



P = w 
Q = uj 



1 



1 



1 1 



\t\ \t-ti\ 



QtxP 



Pt X 



(51) 



where Pi =P^Pi = -P, Q2 ^Q.Qs.^ ~Q- One should 
specify that under construction of approximate solutions 
there are additional basic three vectors by means of which 
solutions at the meeting point of particles are built. One 
should also introduce the appropriate set of times ij^in' 
where formally singular solutions are matched. On a time 
scale t' the solution of the system (|5ip has the form 



Q< - 



cos6' Qt + sm9 {cos[uj{?] — r][)] mp^ 
sin[w(r/ - T][)] np^} , t'^ < t < t" , 
cos 6* Pt + sin0 {cos[aj(C — C2)] ™Qt 
sin[w(C-C2)] nQ^} , t'2<t<t" , 



where 



V 



v - V' , c ^ 

-ln\t-t3\ 



e = - ln|t-t4| , 



(52) 



with the initial data taken on the scale t' . The passage of 
approximate solutions from the scale \t"\ to large positive 
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times is described by the solution of the form 

P> = cos^o Qo + sin 0o {cos[a;(?7 - 77*)] mp^ 

+ sin[w(r/ - 77:)] npj , \t"\ < t < Vt > 

Q> = cos 6*0 Po + sinOo {cos[a;(C - C*)] "^Qo 

- sinHC - Co*)] nq J , 

where the three basic vectors and phase are determined 
by the conditions si = S2 = discussed above. Under 
these conditions one can avoid solutions of the compli- 
cated inverse problem of the restoration of basic triple 
vectors using the initial data on the scale t". For the 
case of two dipoles one can take 

With this choice of the phase the passage to large positive 
times for the first and fourth particle occurs somewhere 
on the scale and on the scale ^3 for the second and 
third particle 



^4 27r 

= ^ _ ' Xo2^—~C ■ (53) 



Accordingly, in the case of two dipoles the passage 
to large positive times is more consistent than in the 
particle-dipole case. Now Eq. p5)) becomes (ci — 1, 
C2 = 1) 



tol 



t3 



2n 

Xol = — 



c = iP"Q") . 



For the particular case considered the system (|44| is re- 
duced to 



Thus the approximate solution on the entire time axis 
has been constructed. 
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